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Abstract 

A method to construct non-dissipative non-Dirac-hermitian relativis- 
tic quantum system that is isospectral with a Dirac-hermitian Hamilto- 
nian is presented. The general technique involves a realization of the 
basic canonical (anti-)commutation relations involving the Dirac matri- 
ces and the bosonic degrees of freedom in terms of non-Dirac-hermitian 
operators, which are hermitian in a Hilbert space that is endowed with a 
pre-determined positive-definite metric. Several examples of exactly solv- 
able non-dissipative non-Dirac-hermitian relativistic quantum systems are 
presented by establishing/employing a connection between Dirac equation 
and supersymmetry. 
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1 Introduction 



An emergent view in the current literature is that a formulation of non-dissipative 
quantum systems is possible without the requirement of the Dirac-hermiticity 1 J 
[21 13]. The quantum systems constructed in this way are indeed hermitian in 
a Hilbert space that is endowed with a positive-definite metric different from 
the identity operator and hence, a modified inner-product. The construction 
of the metric in the Hilbert space is a non-trivial task and is not known for a 
majority of the examples considered in this context. A lack of knowledge of 
the metric makes even an otherwise exactly solvable non-Dirac-hermitian quan- 
tum system incomplete, since no expectation values of physical variables can 
be calculated explicitly. An approach taken in Refs. [4] [5] in this regard was 
to construct non-dissipative non-Dirac-hermitian quantum systems with a pre- 
determined positive-definite metric in the Hilbert space. The quantum systems 
constructed following this approach may serve as prototype examples of non- 
dissipative non-Dirac-hermitian quantum models with a complete description. 
These systems may be used to test different ideas related to the subject (HO [6], 
including validity of any approximate and/or numerical method. 

Recently, non-Dirac-hermitian relativistic quantum systems with entirely 
real spectra have been studied [3 |8j [9j [TO] [IT] fTS] [13]. An optical realization 
of relativistic quantum system with non-Dirac-hermitian interaction has also 
been proposed|13). This proposal assumes a significance within the background 
that ideas emanating from combined parity and time-reversal symmetric optical 
systems have been verified experimentally [HI 115) . A majority of the non-Dirac- 
hermitian relativistic models proposed so far are incomplete in the sense that 
the positive-definite metric and/or the modified inner-product in the Hilbert 
space is not known. 

The purpose of this article is to present a class of non-dissipative non-Dirac- 
hermitian relativistic quantum systems admitting entirely real spectra and uni- 
tary time evolution. The metric and hence, the inner-product in the Hilbert 
space is specified explicitly so as to facilitate the computations of expectation 
values of different physical observables. The general method is an extension 
of the technique employed in Refs. [4] [5] to relativistic systems. In particu- 
lar, a non-Dirac-hermitian realization of Dirac matrices and the bosonic vari- 
ables is used, which are hermitian in a Hilbert space that is endowed with a 
pre-determined positive-definite metric. Apart from the general results, sev- 
eral examples of exactly solvable non-dissipative non-Dirac-hermitian quantum 
systems are presented by employing a connection between supersymmetry and 
Dirac equation. 

2 Preliminaries 

The Pauli matrices a a , position variables x a and the momentum operators p a 
are taken to be hermitian in the Hilbert space Hd with the standard inner- 
product (-|-). A Hilbert-space %ri + that is endowed with the positive-definite 
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metric r; + , 



T} + := exp [—2 J ■ n<f>j , J:=L+^, h-h = l, (</>, n a ) e R, a = 1,2, 3, (1) 

and the inner-product ((•!•))« : = is introduced. The operator L with 

L a = e abc x b p c denotes orbital angular momentum in Hd, where e abc denotes the 
three dimensional Levi-Civita tensor. The operator J can thus be identified as 
the total angular momentum of the system. It may be noted that the metric rj + 
can be decomposed as a direct product of a purely bosonic and purely fermionic 
metrics. In particular, ry+ = r/b ® r\f where rjb := exp(—2L ■ h<fi) and r/f := 
exp(—<7 ■ h(f>). The positivity of the metric follows from the fact that the eigen 
values of the operator J • n are real. With the introduction of the similarity 
operator p, 

p-=Vv+ = ex p{^J- fl( t ) )^ ( 2 ) 

a set of non-Dirac-hermitian matrices E a , position variables X a and the mo- 
mentum operators P a may be introduced as follows [21 HI [5]: 

3 

S a := (p)- 1 a a p = ^R ab a b , 
6=1 
3 

X a := (p)- 1 x a p = Y,R ab x\ 
6=1 

3 

P a : = { P y 1 p a p^^R ab P b , 

6=1 

R ab = n a n b (l-cosh(f>)+d ab cosh(f> + ie abc n c sinh(f>. (3) 

The non-Dirac-hermiticity of E a , X a and P a follows from the fact that (R ab )* ^ 
R ab , where /* denotes the complex conjugate of /. It can be shown by using 
the identity Y? a =i R ab R ac = S bc that E • P = a ■ p. This implies that a ■ p and 
S ■ P are hermitian both in Hd as well as in "H v+ ■ The matrices S a obey the 
same algebra satisfied by the Pauli matrices: 

[E a , E b ] = 2ie abc S c , {E a , E fc } = 2S ab 

S± := 1 (S 1 ± iS 2 ) , {E_,E+} = 1, E 2 ± =0, (4) 

and are hermitian in the vector space . Similarly, the position and the 
momentum operators satisfy the basic canonical commutation relations, 

[X a , X b ] = = [P a , P b ] , [X a 7 P b ] = iS ab , (5) 

and are hermitian in T-L-q + - All other commutators involving E a , X a and P a 
vanish identically. 
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3 Dirac Equation in 1 + 1 dimensions 



In this section, a convenient notation x 1 = x and p 1 = p x will be used to denote 
the position and the momentum operators. A non-Dirac-hermitian relativistic 
Hamiltonian in 1 + 1 dimensions is introduced as follows, 

Hid = cr 2 Px + (tM(.t) + a 3 B(x) + V(x) 
A(x) = M(x)cosh<f) + iP(x)sinh(f>, 

B(x) ee iM (x)sinh(j> - P{x)cosh<p, <p e i?, (6) 

where M(x), P(x) and V(x) are real functions of the coordinate. The complex 
functions A(x) and B(x) can be interpreted as scalar and pseudo-scalar poten- 
tials, respectively. The potential V{x) is the time-component of a Lorentz two- 
vector potential. The space-component of the same two-vector has been taken 
as zero, since it can always be gauged away in 1 + 1 dimensions. The Lorentz co- 
variance of the Dirac equation can be shown by introducing 7 := er 1 , 7 1 := ia 3 . 
The matrix 75 may be defined as 7 s = ia 2 . Various representations of the Dirac 
equation in 1 + 1 dimensions have been used previously [51 HU1 IT31 ITC] . All these 
representations are related to each other through unitary transformations. For 
example, the unitary operator U = exp{— ^er 1 ) relates the representation used 
in this article to the one used in Ref. [13] in the context of optical realization 
of non-Dirac-hermitian relativistic quantum systems. 

The Hamiltonian Hid is non-Dirac-hermitian. With the introduction of the 
hermitian and positive-definite operators, 

PlD :=e-^\ VlD :=p\ D = e-^\ (7) 

the non-Dirac-hermitian Hamiltonian H\d can be mapped to a Dirac-hermitian 
Hamiltonian, 

^id = PidH\dPid 

= a 2 p x + a l M{x) - a 3 P{x) + V{x). (8) 

The similarity operator piD is obtained from Eq. $2$ by choosing L = 0, n 1 = 
= h 3 ,h 2 = 1. The Hamiltonian Hid is hermitian in the Hilbert space Wrim 
that is endowed with the metric r\\D and a modified norm ((-l-))^^ = 
The Hamiltonian H\d and h\D are isospectral. The function M(x) and P{x) 
may be identified as the scalar and the pseudo-scalar interactions, respectively, 
for the Dirac-hermitian Hamiltonian h±D- 

3.1 Case I: Pure scalar interaction 

The eigen value equation for h\D can be solved exactly for specific choices of the 
interactions. The first example considered here is the vanishing vector potential 
and vanishing Pix). The Hamiltonian h\D for this case may be identified as 
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the supercharge of a supersymmetric quantum system with H s :— h 2 D as the 
corresponding supersymmetric Hamiltonian[16 . In particular, 

H s =p 2 x + M 2 - a 3 M', (9) 

where f'(x) denotes the derivative of f(x) with respect to its argument. An 
exhaustive list of M for which H s is exactly solvable is known|16j. The non- 
Dirac hermitian Hamiltonian Hid, 

Hid = ° 2 Vx + a 1 M(x)coshcf) + ia 3 M {x)sinh<f>, (10) 

is thus exactly solvable corresponding to each M for which H s is exactly solvable. 
The energy eigen values E of H\d is determined as E — ±\fW s , where E s 
is the energy of H s . It may be recalled here that E 8 being the energy of a 
supersymmetric theory, it is always semi-positive definite. The eigen-spinor ^ 
of Hid is related to the eigen functions ip of H s as, = Pif, 1 ^- A complete set 
of orthonormal eigen functions ip iiiT-Ld ensures orthonormality of ^ in Urnr, ■ 



3.2 Case II: Scalar & Pseudo-scalar interactions 

The second example describes the case of a vanishing vector potential and non- 
vanishing M(x) and P(x). The square of the operator hiD for this limiting case 
reads, 

h\ D = pi + M 2 + P 2 - a 3 M' - a x P'. (11) 

A unitary transformation U := exp(— i^a 2 ), 9 G [0, 27r] maps h\ D to the Hamil- 
tonian Hi, 



Hi := U- x h 2 1D U 

= pI 

- a 3 (P'sin6 + M'cosO). (12) 



„2 , M 2 + pi + a l ( M > sin Q _ p'^ 



It may be noted that U is unitary 'mH-D as well as in W^ 1D . The term con- 
taining a 1 vanishes identically for a fixed value of 9, when the following relation 
involving the scalar and the pseudo-scalar potential is satisfied: 

P = o,qM + bo, ao,&o S R, 9 = tan~ 1 ao- (13) 

Thus, the Hamiltonian Hi becomes diagonal in this limit and depends on only 
one arbitrary function M(x). 

The Hamiltonian Hi can be re- written in a manifestly supersymmetric form, 

H := Hi — „ 



= pl + W 2 - <t 3 W, (14) 
where the superpotential W(x) is introduced as, 



W(x)^^4M+-^t=. (15) 

\ 1 + <2n 
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Thus, the known relation between supersymmetry and the Dirac equation in 
^dUBEI] is extended to the case of non-Dirac-hermitian quantum systems. 
There is an exhaustive list of superpotentials W(x) for which H is exactly 
solvable|16j. The non-Dirac-hermitian relativistic Hamiltonian, 



Hid = cr 2 Px + cr 1 [e^M+ + e~*M_] + a 3 



M ± ee 1±^-W{x) - (oo T i) , (16) 
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is exactly solvable for each superpotential W(x) for which H is exactly solvable. 
Moreover, a consistent quantum description of Hid, including entirely real spec- 

tra and unitary time-evolution is admissible. An eigen-spinor VP = I ^ j of 

H\d with energy eigenvalues E is related to the corresponding eigen-functions 

" > : ' 



ip=[y)otH with energy eigenvalues E\d through the relation, 
\1p2J 

vf = (p- D U) j,. (17) 

The orthonormality of in Hr) 1D is guaranteed, if the eigenfunctions ip are 
constructed as a complete set of orthonormal eigenfunctions in The eigen- 
values E are determined in terms of the eigenvalues E through the relation, 



E=± f + rh- (18) 

The parameter '60' can be chosen as zero without any loss of generality, since 
E being the energy eigen values of a supersymmctric theory is always semi- 
positive definite. In this limit, E = ±\/e~ and the zero-energy ground state of 
the supersymmetric Hamiltonian H corresponds to the zero modes of the Dirac 
Hamiltonian H id given in Eq. (fTS)) . The choice of a reflection- less superpoten- 
tial W{x) produces a perfectly transparent potential for the non-Dirac-hermitian 
Hamiltonian. 

3.3 Case III: Scalar, Pseudo-scalar &c Vector interactions 

The third example presented in this article is for non- vanishing M(x), P(x) 
and vector potential V(x). The eigen- value equation of hiD with arbitrary 
P(x),M{x) and V(x) reads, 

P(x) + V(x) -i Px + M(x) \ ( rpW \ f \ h q\ 

i Px +M(x) -p(x) + v(x)J \i>(~) ) ~ v^ H y' 

which can be decoupled by introducing the functions 9±(x), x±i x )i 

X±{x) ^mjk' e ± (x) = E±P(x)-V(x), (20) 
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and redefining ip( ±s> in terms of as, 

4>^(x) = e ± -i^(x). (21) 

The decoupled equations have the form of stationary state Schrodinger equa- 
tions, 

H^^E 2 ^, H ± :=~ + V ± (x), 

ax 

V± = M 2 + P 2 -V 2 + 2EV + xl~ X± T (M' - M X ±) ■ (22) 

Unlike the previous two cases, the potentials V± appearing in the effective 
Schrodinger equations explicitly depend on the energy E of the Dirac equation. 
In general, the potentials V± can not be identified as partner superpotentials 
of a supersymmetric Hamiltonian. Such a scenario may arise for very specific 
choices of the interactions, the study of which is beyond the scope of this article. 
Nevertheless, either H + or iJ_ can be embedded in a supersymmetric theory 
and the eigen- value equation can be solved analytically. 

The eigen value equation of the Hamiltonian hiu can be solved by first 
identifying an exactly solvable H- or H + . The eigenfunction of H + can 
then be determined in terms of the eigenfunction ) of if_ or the vice verse, 
by using the equations (p~9|) . (|20|) . and (f2T|) . A choice of the vector potential in 
terms of P{x) as, 

V(x) =P(x) + ai, ai^E, a x € R, (23) 

leads to vanishing \+- With a further choice of P(x) = b\M{x),b\ € R, the 
potential V± reads, 

V+ = W 2 - W +ai(2E- ai ) - 6?(a a - E) 2 , 

V- = W 2 + W + a 1 (2-a 1 )-b 2 1 {a 1 - E) 2 +[x--X-- M X -], 

w . m-m«,-.), (1+2 ^ri:u^ <24) 

Although H + and H- do not become super-partners of each other for this case, 
H+ can be embedded in a supersymmetric theory. A supersymmetric quantum 
system in Hd niay be introduced in terms of the operators A, A' and the 
Hamiltonian H^ 1 \H^> as follows: 

A = ip x + W(x), A 1 = -ip x + W(x) 
d 2 

H {1) := A* A = —r^ + W 2 - W, 
dx z 

HW := Atf = -^- + W 2 + W. (25) 
dx z 

The eigenfunctions and the corresponding eigenvalues of H^> and are 
denoted as ipn\d% an( i ipri^^nK respectively. The Hamiltonian H + and 
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differ by a constant and commute with each other. Thus, the eigen- values E^ 
and the corresponding eigen functions ipn of h\D are determined as, 

E± = ai ±(^^j\ v^v/^F^w. (26) 

The convention to be followed henceforth is that for positive (negative) energy 
solutions E + (E~) will be taken in the expressions of ip+ an d 4>n- ^ follows 
from Eqs. ([19]), (J20]), (|2l|), ([25]) and <(26]) that the eigen functions of h 1D 
corresponding to energy eigen values E^ is determined in terms of ip 1 ^ as, 

Tpt } = r4 [A-hffi-ai)] W 

V E n 
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(27) 



Ol 

where V^ 2 ' = = ^i 2 ! ■ The Hamiltonian .ffW is exactly solvable for several 
choices of W^as) corresponding to shape invariant potentials[l6]. Thus, the 
Hamiltonian h\D and hence, Hid are also exactly solvable for the specific forms 
of V(x), P(x) and M(x) mentioned above for which W(x) produces shape- 
invariant potentials. 

A few comments are in order before the end of this section: 

(i) A choice V(x) = —P(x) + a 2 € R leads to vanishing With a further 
choice of P(x) = b 2 M{x\b2 € R and following the procedure outlined above, 
the complete eigen value equation can be solved analytically by embedding ii_ 
in a supersymmetric theory. 

(ii) The connection between Dirac equation with scalar interaction and super- 
symmetry is well known [TBI 117) . It appears that a connection between super- 
symmetry and Dirac equation with scalar and pseudo-scalar interactions has 
been established for the first time in this article. The method outlined in this 
article for the study of Dirac equation with most general coupling makes an 
indirect connection with supersymmetry for specific reductions of different in- 
teractions. Further investigations along this line may lead to exactly solvable 
H\d with more general interactions. 



4 Dirac Equation in 2 + 1 dimensions 

In 2 + 1 dimensions, a non-dissipative non-Dirac-hermitian relativistic quantum 
system is introduced as follows, 

H 2D = S 1 !! 1 + £ 2 n 2 + <r 3 (j){X 1 ,X 2 ) + ViX 1 .^ 2 ), 

ir a =p a - A a (x 1 ,x 2 ), IF = P a -A a (X\X 2 ), a = 1,2, (28) 

where <p is the scalar potential. The interaction V is the time component of a 
Lorentz vector, while A 1 and A 2 are the space components of the same vector. 
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The similarity operator p2D and the metric r]2D are defined as, 

a 3 

P2D = exp[—<j)(L 3 + -g")]) V2D ■= pl D , (29) 

which can be obtained from Eqs. ([]]) and ^ by choosing n 1 = = n 2 , n 3 = 1. 
A further choice of L 3 = allows to have Rasbha-type spin-orbit interaction 
with purely imaginary coupling in the Hamiltonian Hid- The results of Ref. 
[12] may be reproduced easily in this special limit. The discussion below is for 
the general case of L 3 7^ 0. 

The Hamiltonian H2D is hermitian in the Hilbert space that is endowed 
with the metric T]2D and can be mapped to a Dirac-hermitian Hamiltonian, 

= P2dHp 2 l) 

= o-V +(T 2 7r 2 + cr 3 (j>(x 1 ,x 2 ) + V{x 1 ,x 2 ). (30) 

An analytically solvable h^D corresponds to an exactly solvable H2D- An ex- 
ample in this regard may be the relativistic Landau problem. In particular, 
choosing <\> — 0, A 1 — 0, A 2 — Bx 1 , V — Ex 1 and using the Landau gauge, the 
Hamiltonian h,2D describes a fermion in a crossed uniform electric and magnetic 
fields with magnitude E and B, respectively. Such a Hamiltonian also appears 
in the context of graphene and analytical results are known 18 . 



5 Dirac Equation in 3 + 1 dimensions 

A non-Dirac-hermitian relativistic quantum system in 3 + 1 dimensions may be 
introduced as follows, 

H 3D = a- p + p^(X\X 2 , X 3 ) + i^(X\X 2 , X 3 ) + V{X\X 2 , A 3 ), (31) 

where the Dirac-representation is used for the matrices a, ft and 75: 

The scalar and the pseudo-scalar potentials are denoted by ^(X 1 , X 2 , X 3 ) and 
^(A 1 , A 2 , X 3 ), respectively. The potential term V can be interpreted as the 
time component of a Lorentz four-vector. It may be noted that the non-Dirac- 
hcrmiticity of H^d is due to the presence of A 1 , A 2 , A 3 , which are complex 
functions of their arguments. 

Introducing the similarity operator p^o and the metric operator 773D, 

Psd = (q ° p ^j , V3D ■= pI d , (33) 

a non-Dirac-hermitian realization of the Dirac matrices may be obtained through 
the replacement of a, /3 by cfi, /3i : 

*i:=(g 0)' Pi:=P - (34) 
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The T matrices in this non-Dirac-hermitian realization have the form, 

r :=/?, Ti-Poi, r 5 := 75 . (35) 

It may be noted that a\ ■ P = a ■ p, implying that these two quantities are 
hermitian in Hd as well as in V. V3D . The Hamiltonian can be mapped to a 
Dirac-hermitian Hamiltonian h 3 o as, 

hzD '■= PzdHzdPzd 

= a •p + /3$(x 1 ,x 2 ,x 3 ) + i^(x\x 2 ,x 3 ) + V(x 1 ,x 2 ,x 3 ). (36) 

The Hamiltonian H^£> is hermitian in the Hilbert space H V3D with the modified 
norm ((-\-)) mD - Exactly solvable Dirac-hermitian relativistic quantum systems 
are very few in 3 + 1 dimensions and no attempt has been made in this article 
to present and/or analyze a non-Dirac-hermitian quantum system with such 
properties. 

6 Summary & Conclusions 

A class of non-dissipative non-Dirac-hermitian relativistic quantum systems 
those are isospectral with known Dirac-hermitian Hamiltonian has been con- 
structed in one, two and three space dimensions. A non-Dirac-hermitian re- 
alization of the Dirac matrices and the bosonic operators which are hermitian 
in a Hilbert space that is endowed with a positive-definite metric enables to 
construct this class of relativistic quantum systems. The fundamental canoni- 
cal commutation and anti-commutation relations among different variables re- 
main unchanged for the non-Dirac-hermitian realization. Exact solvability of a 
subclass of these systems has been shown by establishing and/or employing a 
connection between supersymmetry and Dirac equation. An indirect connection 
between Dirac equation with the most general coupling and supersymmetry has 
been established in the process. It is to be seen whether or not any one or more 
non-Dirac-hermitian Dirac equations from the large class of models considered 
in this article have relevance in the realistic physical world. 
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